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Abstract 

We discuss and motivate the form of the generator of a nonlinear 
quantum dynamical group "designed" so as to accomplish a unification 
of quantum mechanics (QM) and thermodynamics. We call this nonrel- 
ativistic theory Quantum Thermodynamics (QT). Its conceptual founda- 
tions differ from those of (von Neumann) quantum statistical mechanics 
(QSM) and (Jaynes) quantum information theory (QIT), but for thermo- 
dynamic equilibrium (TE) states it reduces to the same mathematics, and 
for zero entropy states it reduces to standard unitary QM. By restricting 
the discussion to a strictly isolated system (non-interacting, disentangled 
and uncorrelated) we show how the theory departs from the conventional 
QSM/QIT rationalization of the second law of thermodynamics. The non- 
linear dynamical group of QT is construed so that the second law emerges 
as a theorem of existence and uniqueness of a stable equilibrium state for 
each set of mean values of the energy and the number of constituents. 
To achieve this QT assumes — fc B Trplnp for the physical entropy and is 
designed to implement two fundamental ansatzs. The first is that in addi- 
tion to the standard QM states described by idempotent density operators 
(zero entropy), a strictly isolated system admits also states that must be 
described by non-idempotent density operators (nonzero entropy). The 
second is that for such additional states the law of causal evolution is 
determined by the simultaneous action of a Schrodinger-von Neumann- 
type Hamiltonian generator and a nonlinear dissipative generator which 
conserves the mean values of the energy and the number of constituents, 
and (in forward time) drives the density operator in the 'direction' of 
steepest entropy ascent (maximal entropy increase). The resulting pos- 
itive nonlinear dynamical group (not just a semi-group) is well-defined 
for all nonequilibrium states, no matter how far from TE. Existence and 
uniqueness of solutions of the (Cauchy) initial state problem for all den- 
sity operators, implies that the equation of motion can be solved not only 
in forward time, to describe relaxation towards TE, but also backwards 
in time, to reconstruct the 'ancestral' or primordial lowest entropy state 
or limit cycle from which the system originates. 



1 Introduction 



Several authors have attempted to construct a microscopic theory that includes 
a formulation of the second law of thermodynamics [Tl [21 HI El El El E] - Some 



approaches strive to derive irreversibility from a change of representation of re- 
versible unitary evolution, others from a change from the von Neumann entropy 
functional to other functionals, or from the loss of information in the transi- 
tion from a deterministic system to a probabilistic process, or from the effect of 
coupling with one or more heat baths. 

We discuss the key elements and features of a different non-standard theory 
which introduces de facto an ansatz of "intrinsic entropy and instrinsic irre- 
versibility" at the fundamental level [3110], and an additional ansatz of "steepest 
entropy ascent" which entails an explicit well-behaved dynamical principle and 
the second law of thermodynamics. To present it, we first discuss an essential 
fundamental concept. 

2 States of a strictly isolated individual system 

Let us consider a system A and denote by R the rest of the universe, so that the 
Hilbert space of the universe is Har = Ha®Hr. We restrict our attention to a 
"strictly isolated" system A, by which we mean that at all times, — oo < t < oo, 
A is uncorrelated (and hence disentangled) from R, i.e., par. = Pa®Pr, and 
non-interacting, i.e., Har. — Ha®Ir + Ia^Hr. 

Many would object at this point that with this premise the following discus- 
sion should be dismissed as useless and unnecessary, because no "real" system is 
ever strictly isolated. We reject this argument as counterproductive, misleading 
and irrelevant, for we recall that Physics is a conceptual edifice by which we 
attempt to model and unify our perceptions of the empirical world (physical re- 
ality [H]). Abstract concepts such as that of a strictly isolated system and that 
of a state of an individual system not only are well-defined and conceivable, but 
have been keystones of scientific thinking, indispensable for example to struc- 
ture the principle of causality. In what other framework could we introduce, 
say, the time-dependent Schrodinger equation? 

Because the dominant theme of quantum theory is the necessity to accept 
that the notion of state involves probabilistic concepts in an essential way [12] , 
established practices of experimental science impose that the construct "proba- 
bility" be linked to the relative frequency in an "ensemble" . Thus, the purpose 
of a quantum theory is to regularize purely probabilistic information about the 
measurement results from a "real ensemble" of identically prepared identical 
systems. An important scheme for the classification of ensembles, especially 
emphasized by von Neumann |13j , hinges upon the concept of ensemble "homo- 
geneity" . Given an ensemble it is always possible to conceive of it as subdivided 
into many sub-ensembles. An ensemble is homogeneous iff every conceivable 
subdivision results into sub-ensembles all identical to the original (two sub- 
ensembles are identical iff upon measurement on both of the same physical 
observable at the same time instant, the outcomes yield the same arithmetic 
mean, and this holds for all conceivable physical observables). It follows that 
each individual member system of a homogeneous ensemble has exactly the 
same intrinsic characteristics as any other member, which therefore define the 
"state" of the individual system. In other words, the empirical correspondent 
of the abstract concept of "state of an individual system" is the homogeneous 
ensemble (sometimes also called "pure" [21 [HI [16] or "P r0 P er " [HUH])- 

We restrict our attention to the states of a strictly isolated individual system. 
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By this we rule out from our present discussion all heterogeneous preparations, 
such as those considered in QSM and QIT, which are obtained by statistical 
composition of different homogeneous component preparations. Therefore, we 
concentrate on the intrinsic characteristics of each individual system and their 
irreducible, non-statistical probabilistic nature. 

3 Broader quantum kinematics ansatz 

According to standard QM the states of a strictly isolated individual system are 
in one-to-one correspondence with the one-dimensional orthogonal projection 
operators on the Hilbert space of the system. We denote such projectors by 
the symbol P. If \tjj) is an eigenvector of P such that P\tp) — \ijj) and (4>\4') — 
1 then P = \ip)(ip\. It is well known that differently from classical states, 
quantum states are characterized by irreducible intrinsic probabilities. We need 
not elaborate further on this point. We only recall that — TrPlnP = 0. 

Instead, we adhere to the ansatz [19] that the set of states in which a strictly 
isolated individual system may be found is broader than conceived in QM, specif- 
ically that it is in one-to-one correspondence with the set of linear operators p 
on H, with p^ = p, p > 0, Trp = 1, without the restriction p 2 = p. We call these 
the "state operators" to emphasize that they play the same role that in QM is 
played by the projectors P, and that they are associated with the homogeneous 
preparation schemes. This fundamental ansatz has been first proposed by Hat- 
sopoulos and Gyftopoulos [T!5] . It allows an implementation of the second law of 
thermodynamics at the fundamental level in which the physical entropy, given 
by s(p) — — fc B Trplnp, emerges as an intrinsic microscopic and non-statistical 
property of matter, in the same sense as the (mean) energy e(p) = TrpH is an 
intrinsic property. 

We first assume that our isolated system is an indivisible constituent of 
matter, i.e., one of the following: 

• A single strictly isolated d-level particle, in which case TC = Hd = ®k=o^ e k 
where et is the A:-th eigenvalue of the (one-particle) Hamiltonian Hi and 
TL £k the corresponding eigenspace). Even if the system is isolated, we do 
not rule out fluctuations in energy measurement results and hence we do 
not assume a "microcanonical" Hamiltonian (i.e., H = e^Pft^ f° r some 

k) but we assume a full "canonical" Hamiltonian H = Hi = J^k efc -^e fe • 

• A strictly isolated ideal Boltzmann gas of non-interacting identical in- 
distinguishable d-level particles, in which case H is a Fock space, TL = 
3~d = ©^o^f Again, we do not rule out fluctuations in energy nor in 
the number of particles, and hence we do not assume a canonical num- 
ber operator (i.e., N = zP H ®z for some z) but we assume a full grand 

d 

canonical number operator ./V = Y] n —n TiR H ® n and a full Hamiltonian 

d 

H = Y^=v H n P n »n where H n = Y^j=i{H^)j®^~f 1S the n-particlc Hamil- 

d ' 

tonian on Hf n , {H\)j denotes the one-particle Hamiltonian on the J-th 
particle space {Hd).j and I-j the identity operator on the direct product 
space <8>J-_i k^jO^<i)k of all other particles. Note that [H, N] = 0. 

• A strictly isolated ideal Fermi-Dirac or Bose-Einstein gas of non-interacting 
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identical indistinguishable d-level particles, in which case ri is the anti- 
symmetric or symmetric subspace, respectively, of the Boltzmann Fock 
space just defined. 

We further fix ideas by considering the simplest quantum system, a 2-level 
particle, a qubit. It is well known |20j that using the 3-vector cr — (ai, ui, 03) 
of Pauli spin operators, [o"j,<Tfe] = tju^i, we can represent the Hamiltonian 
operator as H — huo (§/ + h • cr) where h is a unit-norm 3-vector of real scalars 
(hi, fi2, /13), and the density operators as p = \I + r • cr where r is a 3-vector 
of real scalars (r\,ri,r^) with norm r = |r| < 1, and r = 1 iff p is idempotent, 
p 2 = p. 

If the 2-level particle is strictly isolated, its states in standard QM are one- 
to-one with the unit-norm vectors ip in Ti or, equivalently, the unit-trace one- 
dimensional projection operators on Ti, P^ = {^(iplip) -1 i.e., the idempo- 
tent density operators p 2 = p. Hence, in the 3-dimensional euclidean space 
(f\, T2, Tz), states map one-to-one with points on the unit radius 2-dimensional 
spherical surface, r = 1, the "Bloch sphere". The mean value of the energy is 
e(p) = TrpH = i(l + h • r) and is clearly bounded by < e(p) < %u. The 
set of states that share a given mean value of the energy are represented by the 
1-dimcnsional circular intersection between the Bloch sphere and the constant 
mean energy plane orthogonal to h defined by the h • r = const condition. The 
time evolution according to the Schrodinger equation ip = —iHip/h or, equiv- 
alently, P^ — —i[H,P^]/h or [20] r = wh x r yields a periodic precession of 
r around h along such 1-dimensional circular path on the surface of the Bloch 
sphere. At the end of every (Poincare) cycle the strictly isolated system passes 
again through its initial state: a clear pictorial manifestation of the reversibility 
of Hamiltonian dynamics. 

At the level of a strictly isolated qubit, the Hatsopoulos-Gyftopoulos ansatz 
amounts to accepting that the two-level system admits also states that must be 
described by points inside the Bloch sphere, not just on its surface, even if the 
qubit is noninteracting and uncorrelated. The eigenvalues of p arc (1 ± r)/2, 
therefore the isoentropic surfaces are concentric spheres, 

/x , /l + r . 1 + r 1-r, l-r\ 
s(p) = s(r) = -k B f — In + In —\ . (1) 

The highest entropy state with given mean energy is at the center of the disk ob- 
tained by intersecting the Bloch sphere with the corresponding constant energy 
plane. Such states all lie on the diameter along the direction of the Hamiltonian 
vector h and are thermodynamic equilibrium (maximum entropy principle [21j). 

Next, we construct our extension of the Schrodinger equation of motion valid 
inside the Bloch sphere. By assuming such law of causal evolution, the second 
law will emerge as a theorem of the dynamics. 

4 Steepest-entropy-ascent ansatz 

Let us return to the general formalism for a strictly isolated system. We go 
back to the qubit example at the end of the section. 

As a first step to force positivity and hermiticity of the state operator p we 
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assume an equation of motion of the form 

^= P E(p) + tf(p)p = ^-p(^5E(p)) + (y/pE{p)fjp, (2) 

where E(p) is a (non-hermitian) operator- valued (nonlinear) function of p that 
we call the "evolution" operator. Without loss of generality, we write E — 
E + + iE_ where E + = (E + £t)/2 and E- = (E — E^)/2i are hermitian 
operators, so that Eq. @ takes the form 

^ [ = - l [E_(p),p]+{E + {p),p} , (3) 

with [ • , • ] and { • , • } the usual commutator and anti-commutator, respectively. 

We consider the space of linear (not necessarily hermitian) operators on Ti 
equipped with the real scalar product 

(F\G) = Tr(F*G + G f F)/2 , (4) 

so that for any time-independent hermitian observable R on Ti., the rate of 
change of the mean value r{p) = Tr(pR) = (^p\^/pR) can be written as 

^ = Tr(^i?) = 2 UpE\ ^pR) , (5) 

from which it follows that a set of Ti (p) 's is time invariant iff ^fp~E is orthogonal 
to the linear span of the set of operators *J~pRi, that we denote by C{^fpRi\. 

For an isolated system, we therefore require that, for every p, operator ^fp~E 
be orthogonal [in the sense of scalar product (|H)] to the linear manifold C{^fp~Ri} 
where the set {y/pRi} always includes yfpl, to preserve Trp = 1, and t/p~H, to 
conserve the mean energy e(p) = TrpH. For a field of indistinguishable particles 
we also include y/pN to conserve the mean number of particles n(p) = Tr pN . 
For a free particle we would include yfpP x , \fpPyi \fpP% to conserve the mean 
momentum vector p(p) = TrpP, but here we omit this case for simplicity [50] . 

Similarly, the rate of change of the entropy functional can be written as 

^- = (^pE\-2k B [^p + y /phxp]) , (6) 

where the operator —2k B [y^+y^lnp] may be interpreted as the gradient 
(in the sense of the functional derivative) of the entropy functional s(p) — 
— k B Trphip with respect to operator ^fp (for the reasons why in our theory 
the physical entropy is represented by the von Neumann functional, see Refs. 

mm)- 

It is noteworthy that the Hamiltonian evolution operator 

E H = iH/h , (7) 

is such that y/pE H is orthogonal to C{y/pl, y/pH(, ^/pN)} as well as to the en- 
tropy gradient operator —2k B [y/p + ^/p\n p\ . It yields a Schrodinger-Liouville- 
von Neumann unitary dynamics 

^ = pE H + El I p = -^[H,p} , (8) 
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which maintains time-invariant all the eigenvalues of p. Because of this feature, 
all time-invariant (equilibrium) density operators according to Eq. ||5J) (those 
that commute with H) are globally stable [53] with respect to perturbations 
that do not alter the mean energy (and the mean number of particles). As a 
result, for given values of the mean energy e(p) and the mean number of particles 
n(p) such a dynamics would in general imply many stable equilibrium states, 
contrary to the second law requirement that there must be only one (this is the 
well-known Hatsopoulos-Keenan statement of the second law [21] , which entails 
[2"T] the other well-known statements by Clausius, Kelvin, and Caratheodory). 

Therefore, we assume that in addition to the Hamiltonian term E H , the 
evolution operator E has an additional component E 



D ' 



E 



E T 



(9) 



that we will take so that y/pE D is at any p orthogonal both to y/pE H and 
to the intersection of the linear manifold C{^fpRi] with the isoentropic hyper- 
surface to which p belongs (for a two level system, such intersection is a one- 
dimensional planar circle inside the Bloch sphere). In other words, we assume 
that ^/~p~E D is proportional to the component of the entropy gradient operator 
— 2k B [y/p + y^lnp] orthogonal to C{^/pRi}, 



1 



D 



2t( P ) 



(10) 



where we denote the "constant" of proportionality by l/2r(/?) and use the fact 
that yfp has no component orthogonal to £{^/pI, y/pH(, ^fp~N)}. 

It is important to note that the "intrinsic dissipation" or "intrinsic relax- 
ation" characteristic time r(p) is left unspecified in our construction and need 
not be a constant. All our results hold as well if r(p) is some reasonably well 
behaved positive definite functional of p. The empirical and/or theoretical de- 
termination of r(p) is a most challenging open problem in our research program. 
For example, it has been suggested [25] that the experiments by Franzen [26] 
(intended to evaluate the spin relaxation time constant of vapor under vanishing 
pressure conditions) and by Kukolich [37J (intended to provide a laboratory vali- 
dation of the time-dependent Schrodinger equation) both suggest some evidence 
of an intrinsic relaxation time. 

Using standard geometrical notions, we can show [9] [TO] EH [29] that given 
any set of linearly independent operators {^/pRi} spanning L{^fpl, ^/pH(, t/p~N)} 
the dissipative evolution operator takes the explicit expression 
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-VpAM(p) 



D 2k B T(p) 

where M(p) is a "Massieu- function" operator defined by 



(11) 



M(p) - 
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(12) 



and we use the notation (F and G hermitian) 



S = -k B P Ranp \np , (13) 
AF — F ~ Tt(pF)I , (14) 

(AF AG) = (VpAF|VpAG)-iTr(p{AF,AG}) , (15) 
r({«/pRi}) = det[{ARiARj)} (a Gram determinant) . (16) 

The Massieu-function operator defined by Eq. (|12p generalizes to any non- 
equilibrium state the well-known equilibrium Massieu characteristic function 
s(pte) — P z(pte) [+/3pn(pTE)]' As a result, our full equation of motion takes 
the form 

Ghcorghiu-Svirschevski |30j re-derived our nonlinear equation of motion 
from a variational principle that in our notation may be cast as follows |29] , 

max su bject to = and (JpE D \JpE D ) = c 2 {p) , (18) 

Vp e d ™ ™ 

where ro(p) = Trp, ri(p) = Tr Hp [, ^(p) = TriVp], and c 2 (p) is some positive 
functional. The last constraint means that we are not really searching for "max- 
imal entropy production" but only for the direction of steepest entropy ascent, 
leaving unspecified the rate at which such direction "attractsd" the state of the 
system. The necessary condition in terms of Lagrange multipliers is 

and, using Eqs. §5$ and © becomes 

-2k*{ y fp+^p\ixp)-2^\ i yfpR i -2\ Q yfpE D =Q , (20) 

which inserted in the constraints and solved for the multipliers yields Eq. pip . 

The resulting rate of entropy change (entropy generation by irreversibility, 
for the system is isolated) is given by the equivalent expressions 

MP) = -kJ^^=Ak B r(p)( V pE D \^pE D ) (21) 



dt dt 

k B r(051np,{ypfl t }) 1 T(^pS,{VpR t }) 



> . (22) 



Because a Gram determinant T(y/pXi, . . . , y/pXjy) = det[(AX, AX,)] is either 
strictly positive or zero iff operators {^Jp~Xi\ are linearly dependent, the rate 
of entropy generation is either a positive semi-definite nonlinear functional of p, 
or it is zero iff operators ^fpS, y/pl, y/pH(, ,J~pN) are linearly dependent, i.e., iff 
the state operator is of the form 

BeM-PH(+ v N)]B 
P Tr{Bexp[-pH(+vN)]) ' 1 j 
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for some "binary" projection operator B (B 2 = B, eigenvalues either or 1) 
and some real scalar (s) j3 (and v). Non-dissipative states are therefore all and 
only the density operators that have the nonzero eigenvalues "canonically" (or 
"grand canonically") distributed. For them, y/pE D = and our equation of 
motion (|17|) reduces to the Schrodinger-von Neumann form ihp = [H, p]. Such 
states are either equilibrium states, if [B, H) = 0, or belong to a limit cycle and 
undergo a unitary hamiltonian dynamics, if [B, H] =/= 0, in which case 

p{t) = B{t) exp[-f3H (+vN)]B{t)/Tr[B{t) exp[-0H (+vN)}\ , (24) 
B{t) = U{t)B{Q)U- l (t) , U{t) = exp{-itH/h) . (25) 

For TrB = 1 the states (f2"3")l reduce to the (zero entropy) states of standard 
QM, and obey the standard unitary dynamics generated by the usual time- 
dependent Schrodinger equation. For B = I we have the maximal-entropy 
(thermodynamic-equilibrium) states, which turn out to be the only globally sta- 
ble equilibrium states of our dynamics, so that the Hatsopoulos-Keenan state- 
ment of the second law emerges as an exact and general dynamical theorem. 

Indeed, in the framework of our extended theory, all equilibrium states and 
limit cycles that have at least one null eigenvalue of p are unstable. This is be- 
cause any neighboring state operator with one of the null eigenvalues perturbed 
(i.e., slightly "populated") to a small value e (while some other eigenvalues are 
slightly changed so as to ensure that the perturbation preserves the mean energy 
and the mean number of constituents), would eventually proceed "far away" to- 
wards a new partially-maximal-entropy state or limit cycle with a canonical 
distribution which fully involves also the newly "populated" eigenvalue while 
the other null eigenvalues remain zero. 

It is clear that the canonical (grand-canonical) density operators pte = 
exp[— j3H (+^7V)]/Tr(exp[— f}H (+vN)]) are the only stable equilibrium states, 
i.e., the TE states of the strictly isolated system. They are mathematically 
identical to the density operators which also in QSM and QIT are associated 
with TE, on the basis of their maximizing the von Neumann indicator of sta- 
tistical uncertainty — Trplnp subject to given values of TrHp (and TrNp). Be- 
cause maximal entropy mathematics in QSM and QIT successfully represents 
TE physical reality, our theory, by entailing the same mathematics for the stable 
equilibrium states, preserves all the successful results of equilibrium QSM and 
QIT. However, within QT such mathematics takes up an entirely different phys- 
ical meaning. Indeed, each density operator here does not represent statistics 
of measurement results from a "heterogeneous" ensemble, as in QSM and QIT 
where, according to von Neumann's recipe |13U31j . the "intrinsic" uncertainties 
(irreducibly introduced by standard QM) are mixed with the "extrinsic" un- 
certainties (related to the heterogeneity of its preparation, i.e., to not knowing 
the exact state of each individual system in the ensemble). In QT, instead, 
each density operator, including the maximal-entropy stable TE ones, repre- 
sents "intrinsic" uncertainties only, because it is associated with a homogeneous 
preparation and, therefore, it represents the state of each and every individual 
system of the homogeneous ensemble. 

We noted elsewhere [33] that the fact that our nonlinear equation of motion 
preserves the null eigenvalues of p, i.e., conserves the cardinality dimKer(p) 
of the set of zero eigenvalues, is an important physical feature consistent with 
recent experimental tests (see the discussion of this point in Ref. [30] and ref- 
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erences therein) that rule out, for pure (zero entropy) states, deviations from 
linear and unitary dynamics and confirm that initially unoccupied eigenstates 
cannot spontaneously become occupied. This fact, however, adds nontrivial ex- 
perimental and conceptual difficulties to the problem of designing fundamental 
tests capable, for example, of ascertaining whether decoherence originates from 
uncontrolled interactions with the environment due to the practical impossibil- 
ity of obtaining strict isolation, or else it is a more fundamental intrinsic feature 
of microscopic dynamics requiring an extension of QM like the one we propose. 

For a confined, strictly isolated d-level system, our equation of motion for 
non-zero entropy states (p 2 ^ p) takes the following forms [20l [34]. If the 
Hamiltonian is fully degenerate [H = el, e(p) = e for every p], 



-£ = -^[H,p]--(p\np-pTrplnp) , 
at n t 



(26) 



while if the Hamiltonian is nondegenerate, 

p\np 



d P i \u l 1 



Trp In p 



\{H,P\ 
TrpH 



TrpH In p TrpH TrpH 2 



TrpH 2 - {TrpH) 2 



(27) 



In particular, for a non-degenerate two-level system, it may be expressed in 
terms of the Bloch sphere representation (for < r < 1) as [30] 



ujh x r — — 

T 



1 — r 2 1 — r 
In 



2r 



1 



hxrxh 
1- (h-r) 2 



(28) 



from which it is clear that the dissipative term lies in the constant mean energy 
plane and is directed towards the axis of the Bloch sphere identified by the 
Hamiltonian vector h. The nonlinearity of the equation does not allow a general 
explicit solution, but on the central constant-energy plane, i.e., for initial states 
with r • h = 0, the equation implies [2D] 



d , 1 — r 1 1 — r 

— In = In 

dt 1 + r t 1 + r 



(29) 



which, if r is constant, has the solution 



r(t) = tanh 



, t\ i - Kg) 



(30) 



This, superposed with the precession around the hamiltonian vector, results in a 
spiraling approach to the maximal entropy state (with entropy k B In 2). Notice, 
that the spiraling trajectory is well-defined and within the Bloch sphere for all 
times — oo < t < +oo, and if we follow it backwards in time it approaches as 
t — > — oo the limit cycle which represents the standard QM (zero entropy) states 
evolving according to the Schrodinger equation. 

This example shows quite explicitly a general feature of our nonlinear equa- 
tion of motion which follows from the existence and uniqueness of its solutions 
for any initial density operator both in forward and backward time. This feature 
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is a consequence of two facts: (1) that zero eigenvalues of p remain zero and 
therefore no eigenvalue can cross zero and become negative, and (2) that Trp is 
preserved and therefore if initially one it remains one. Thus, the eigenvalues of p 
remain positive and less than unity. On the conceptual side, it is also clear that 
our theory implements a strong causality principle by which all future as well 
as all past states are fully determined by the present state of the isolated sys- 
tem, and yet the dynamics is physically (thermodynamically) irreversible. Said 
differently, if we formally represent the general solution of the Cauchy problem 
by pit) = A t/ o(0) the nonlinear map A t is a group, i.e., A t+U — A t A u for all t 
and it, positive and negative. The map is therefore "invertible" , in the sense 
that A_ t = A^ 1 , where the inverse map is defined by p(0) = A^ 1 p(t). 

It is a nontrivial observation that the non-invertibility of the dynamical map 
is not at all necessary to represent a physically irreversible dynamics. Yet, 
innumerable attempts to build irreversible theories start from the assertion that 
in order to represent thermodynamic irreversibility the dynamical map should 
be non-invertible. The arrow of time in our view is not to be sought for in 
the impossibility to retrace past history, but in the spontaneous tendency of 
any physical system to internally redistribute its energy (and, depending on the 
system, its other conserved properties such number of particles, momentum, 
angular momentum) along the path of steepest entropy ascent. 



5 Onsager reciprocity 

The intrinsically irreversible dynamics entailed by the dissipative (non-hamiltonian) 
part of our nonlinear equation of motion also entails an Onsager reciprocity the- 
orem. To see this, we first note that any density operator p can be written as 

ED 



/?cx|h V / f .V f i/? 

P TvBeM-EjfjXj) ' 1 j 

where the possibly time-dependent Boolean B is such that B = -Pa an p (= I — 
-fkerp) and the time-independent operators Xj together with the identity I form 
a set such that their restrictions to H' = BH, {I',X'A span the real space of 
hermitian operators on H' = BH. Hence, 

V^lnp = -/oVP ~ Ej fiVpXj > (32) 

x j (p)=Tr(pX j ) , (33) 

s(p) = k B f + k B J2j fj Xj (p) , (34) 
ds(p) 



where k B fj 



(35) 



dxj{p) 

may be interpreted as a "generalized affinity" or force. Defining 

5^=2(2^^). 06) 
as "the dissipative rate of change" of the mean value Xj (p) , we find 



Dt 
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where the coefficients Lij (p) (nonlinear in p) may be interpreted as "generalized 
conductivities" and are given explicitly (no matter how far p is from TE) by 



Lij{p) = Lji(p) 



r(p) 



([VP X i}±C{^/pR t }\ [VP X j]±C{y/pIU}) ( 38 ) 





(AXiAXj) 


(AfliAXj) ■ • • 


(AR k AXj) ■ ■ ■ 




(AXiARi) 


(AfliAHi) ■ • • 


(AR k AR!) ■ ■ ■ 


1 


(AXiAR k ) 


(ARxARt,) ■ ■ ■ 


(AR k AR k ) ■ ■ ■ 


r(p) 









(39) 



and therefore form a symmetric, non-negative definite Gram matrix [Lij(p)], 
which is strictly positive iff all operators [^/pXi]j_c{^/pRi} are linearly indepen- 
dent. 

The rate of entropy generation may be rewritten as a quadratic form of the 
generalized affinities, 



ds(p) 
dt 



(40) 



If all operators [^/pXi]±c{^/pRi} are linearly independent, det[Ly(p)] ^ and 
Eq. (f37|) may be solved to yield 



Dt 



(41) 



and the rate of entropy generation can be written also as a quadratic form of 
the dissipative rates 



dt 



« 3 



Dt Dt 



(42) 



6 Composite systems and reduced dynamics 

The composition of the system is embedded in the structure of the Hilbert space 
as a direct product of the subspaces associated with the individual elementary 
constituent subsystems, as well as in the form of the Hamiltonian operator. In 
this section, we consider a system composed of distinguishable and indivisible 
elementary constituent subsystems. For example: 

• A strictly isolated composite of r distinguishable d-level particles, in which 
case H = ® r J= {Hdj and H = Y?j=i(Hi)j Ij + V where V is some 
interaction operator over H. 

• A strictly isolated ideal mixture of r types of Boltzmann, Fermi-Dirac or 
Bose-Einstein gases of non-interacting identical indistinguishable <ij-level 
particles, J = 1, . . . ,r, in which case Tt is a composite of Fock spaces 
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W = ® r J= iFdj = ®Z=Q ■ ■ ■ ®n r =0 ^fr 1 ® ■ ■ * ® W| ^ where the fact0r 

Fock spaces belonging to Fermi-Dirac (Bose-Einstein) components are re- 
stricted to their antisymmetric (symmetric) subspaces. Again, we assume 
full grand-canonical number operators Nj — X^ 7 =o n jP-H^ n J an d Hamil- 

tonian H = Y,j=i ^7® 2^=0 H njP H ®nj + V. 

d .i 

For compactness of notation we denote the subsystem Hilbert spaces as 

H = H^H 2 ® ■ ■ ■ ®H r = H J ®H 7 , (43) 

where J denotes all subsystems except the J-th one. The overall system is 
strictly isolated in the sense already defined, and the Hamiltonian operator 

r 

H = Y. H .J® I 1+ V > ( 44 ) 
J=i 

where H j is the Hamiltonian on 7i associated with the J-th subsystem when 
isolated and V (on H.) the interaction Hamiltonian among the r subsystems. 

The subdivision into elementary constituents, considered as indivisible, and 
reflected by the structure of the Hilbert space Tt as a direct product of subspaces, 
is particularly important because it defines the level of description of the sys- 
tem and specifies its elementary structure. The system's internal structure we 
just defined determines the form of the nonlinear dynamical law proposed by 
this author [2"B"1 US EI] to implement the steepest entropy ascent ansatz in a 
way compatible with the obvious self-consistency "separability" and "locality" 
requirements [33j . It is important to note that, because our dynamical principle 
is nonlinear in the density operator, we cannot expect the form of the equation 
of motion to be independent of the system's internal structure. 

The equation of motion that we "designed" in [23 EZ] so as to guarantee all 
the necessary features (that we list in Ref. EH]), is 

t = ^^t^^MY^^ , (45) 

where we use the notation [see Ref. [29] for interpretation of (S) J and (H) J ) 
Vp]( AM Ap)) J = [VP](Sy] (46) 

C iVP]( R t j) J } = ^■^oi^fT ] I J ,^fT ] (H j y\^p- J {N kJ ) J (47) 

(FjlGJj = Tr j(F^Gj + G\Fj)/2, (48) 

(R u y' = Tt-KI^RJ , (49) 

(S) J = Tt-KIj^S] , (50) 

and the "internal redistribution characteristic times" r J (p) , s are some positive 
constants or positive functionals of the overall system's density operator p. 

All the results found for the single constituent extend in a natural way to 
the composite system. For example, the rate of entropy change becomes 

d S (p) i r(^j(sy, {^{R u y}) 

at ^fer» T({^pj(R u y}) 



12 



The dynamics reduces to the Schrodinger-von Neumann unitary Hamiltonian 
dynamics when, for each J, there are multipliers X iJ such that 

i 

The equivalent variational formulation is 

/ m p x x sub J ect t0 ~!f^ = and (VPj E dj\VPj E dj)j = c2 j(p) > 

(53) 

where ro(p) — Trp, r\{p) = TrHp [, ^(p) = TtNp], and c^(p) are some posi- 
tive functionals of p. The last constraints, one for each subsystem, mean that 
each subsystem contributes to the overall evolution (for the dissipative non- 
hamiltonian part) by pointing towards its "local perception" of the direction 
of steepest (overall) entropy ascent, each with an unspecified intensity (which 
depends on the values of the functionals Cj (p) , that are inversely related to the 
internal redistribution characteristic times Tj(p)). 

If two subsystems A and B are non-interacting but in correlated states, the 
reduced state operators obey the equations 



dpA 
dt 



l\H A ,p A ] + i- p ^— ^-y{(AM J (p)) J , Pj}®{pa)j , (54) 



JEA 



IT - -^PB^ + it^i^MA^^PjMPBh, (55) 



j- 

JeB 



where (pa)j = Tr 7 (p^), {pb)j — ^j(pb), and operators (AMj(p)) J result 
independent of Hb for every J € A and independent of Ha for every J £ 
B. Therefore, all functionals of pa (local observables) remain unaffected by 
whatever change in B 7 i.e., locality problems are excluded. 



7 Concluding remarks 

According to QSM and QIT, the uncertainties that are measured by the physical 
entropy, are to be regarded as either extrinsic features of the heterogeneity 
of an ensemble or as witnesses of correlations with other systems. Instead, 
we discuss an alternative theory, QT, based on the Hatsopoulos-Gyftopoulos 
fundamental ansatz [19l [31] that also such uncertainties are irreducible (and 
hence, "physically real" and "objective" like standard QM uncertainties) in 
that they belong to the state of the individual system, even if uncorrelated and 
even if a member of a homogeneous ensemble. 

According to QT, second law limitations emerge as manifestations of such 
additional physical and irreducible uncertainties. The Hatsopoulos-Gyftopoulos 
ansatz not only makes a unified theory of QM and Thermodynamics possible, 
but gives also a framework for a resolution of the century old "irreversibility 
paradox", as well as of the conceptual paradox [3T] about the QSM/QIT inter- 
pretation of density operators, which has preoccupied scientists and philosophers 
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since when Schrodinger surfaced it in Ref . [32] . This fundamental ansatz seems 
to respond to Schrodinger prescient conclusion in Ref. [32] : ". . .in a domain 
which the present theory (Quantum Mechanics) does not cover, there is room 
for new assumptions without necessarily contradicting the theory in that region 
where it is backed by experiment." 

QT has been described as "an adventurous scheme" [3B], and indeed it re- 
quires quite a few conceptual and interpretational jumps, but (1) it does not 
contradict any of the mathematics of either standard QM or TE QSM/QIT, 
which are both contained as extreme cases of the unified theory, and (2) for 
nonequilibrium states, no matter how "far" from TE, it offers the structured, 
nonlinear equation of motion proposed by this author which models, determinis- 
tically, irreversibility, relaxation and decoherence, and is based on the additional 
ansatz of steepest-entropy-ascent microscopic dynamics. 

Many authors, in a variety of contexts |35j . have observed in recent years 
that irreversible natural phenomena at all levels of description seem to obey a 
principle of general and unifying validity. It has been named [35] "maximum 
entropy production principle", but we note in this paper that, at least at the 
quantum level, the weaker concept of "attraction towards the direction of steep- 
est entropy ascent" EH [25] is sufficient to capture precisely the essence of the 
second law. 

We finally emphasize that the steepest-entropy-ascent, nonlinear law of mo- 
tion we propose, and the dynamical group it generates (not just a semi-group), 
is a potentially powerful modeling tool that should find immediate application 
also outside of QT, namely, regardless of the dispute about the validity of the 
Hatsopoulos-Gyftopoulos ansatz on which QT hinges. Indeed, in view of its 
well-defined and well-behaved general mathematical features and solutions, our 
equation of motion may be used in phenomenological kinetic and dynamical the- 
ories where there is a need to guarantee full compatibility with the principle of 
entropy non-decrease and the second-law requirement of existence and unique- 
ness of stable equilibrium states (for each set of values of the mean energy, of 
boundary-condition parameters, and of the mean amount of constituents). 
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